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In this paper we report a theoretical study of the static 
dipole polarizability of two one-dimensional structures: (a) 
linear carbon chains C„(n = 2 — 10) and (b) ladder- like 
planar boron chains B n (n = 4 — 14). The polarizabil- 
ities of these chains are calculated both at the Hartree- 
Fock and the correlated level by applying accurate ab ini- 
tio quantum chemical methods. Methods such as restricted 
Hartree-Fock, multi-configuration self-consistent field, multi- 
reference configuration-interaction method, M0ller-Plesset 
second-order perturbation theory, and coupled-cluster singles, 
doubles and triples level of theory were employed. Results ob- 
tained from ab initio wave-function-based methods are com- 
pared with the ones obtained from the density-functional the- 
ory. For the clusters studied, directionally averaged polariz- 
ability per atom for both the systems is seen to increase with 
the chain size. 



I. INTRODUCTION 

The role that atomic and molecular clusters will play 
in future nanotechnologies is indisputable.Ela The ex- 
perimental progress in this field has been breathtaking, 
and novel applications have been found in, areas such 
as molecular transport and optoelectronics.™™ However, 
theoretical research in this area can also play a very 
important role in that, by undertaking calculations on 
clusters of different types, it can help the experimen- 
talists in identifying novel structures for investigation. 
Ab initio calculations on structural and electronic prop- 
erties of atomic clusters are frequently performed, and 
the results are put to test in the experiments. □ One such 
property of clusters is their electric-dipole polarizability 
whose experimental and .theoretical determination is an 
area of intense research.EI The measurements of dipole 
polarizability are frequently used by the experimental- 
ists to characterize the nature of atomic and molecu- 
lar species.EJ It describes the response of the electron 
cloud of the given molecular system to the presence of 
a d.c. electric field, and thus is easily amenable to ex- 
periments. Since the static polarizabilities are the zero- 
frequency limits of the corresponding dynamic quanti- 
ties, they also provide information about the response 
of the system to off-resonant a.c. fields. Most of the 



theoretical calculations of both the structural and elec- 
tronic properties such as static polarizabilities of atomic 
clusters are performed within the framework of the den- 
sity functional theory (DFT). Despite the fact that DFT 
has enjoyed indisputable success in solid state physics 
and quantum chemistry as a computationally cheap rou- 
tine tool for large-scale investigations, it has the draw- 
back that results depend highly on the chosen functional, 
and cannot be improved in a systematic way. Wave- 
function-based ab initio quantum-chemical techniques 
on the other hand are free from this flaw, and provide 
a large array of methods of different accuracy and com- 
putational cost. Moreover, the prediction of reliable val- 
ues of dipole polarizabilities and hyperpolarizabilities by 
rigorous quantum chemical methods has made signifi- 
cant contributions, and added new vigor, to the search 
of novel optical materialsM. Thus, in order to obtain 
reliable estimates for dipole polarizabilities, and also to 
cross check the DFT-based results, it is worthwhile to in- 
vestigate the electron correlation effects in a systematic 
way by using the quantum-chemical many-body tech- 
niques. In this work, we present fully size-consistent 
ab initio calculations to the static dipole polarizability 
of linear carbon clusters C n (n = 2 — 10) and chain- 
like boron clusters B2„(?i = 2 — 7), of increasing size. 
The reason behind our focus on one-dimensional struc- 
tures is that quantum confinement due to reduced di- 
mensions, combined with the possible delocalization of 
the electrons along the backbone, can lead to enhanced 
linear and nonlinear susceptibilities of these structures, 
as compared to their three-dimensional counterparts. In 
the present study, the electron-correlation effects have 
been taken into account by various size-consistent meth- 
ods: multi-reference configuration interaction, second- 
order M0ller-Plesset perturbation theory, coupled-cluster 
singles and doubles (CCSD), and coupled-cluster sin- 
gles and doubles with the perturbative treatment of the 
triples (CCSD(T)). All earlier calculations on these sys- 
tems, with the exception of B4, were performed within 
the framework of DFT, with which we compare our re- 
sults. Next, we briefly review the state-of-the-art of re- 
search on these two types of clusters. 

Carbon clusters have been the subject of research for 
decades as possible key materials for future nanotech- 
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nologiescl. For the smaller systems, up to and including 
those containing nine atoms, linear neutral, positively 
and negatively charged clusters are generated and de- 
tected in experimental. The structures and energetics 
of linear carbon clusters ar£._well studied by employing 
coupled-cluster approaches. Era Recently, we studied the 
ground state of an infinite carbon chain at the ab ini- 
tio level using various many-body approaches including 
CCSD(T).EI Lou et al. have investigated the influence of 
an electricJield on the energetic stability of linear car- 
bon chainsE!! Recently FuentealbaO has calculated the 
static dipole polarizabilities of carbon chain using density 
functionals of the hybrid type in combination with the 
finite-field method. He showed that dipole polarizabil- 
ities are an important quantity for the identification of 
clusters with different numbers of atoms and even for the 
separation of isomers. Here we compare our many-body- 
methods baaed polarizabilities with those computed by 
FuentealbaO using the DFT approaches. 

Boron is a trivalent element with the valence shell 
configuration s 2 p 1 . Although compared to carbon, the 
valence shell of boron is electron deficient, it still ex- 
hibits sp 2 -hybridization with strong directional chemi- 
cal bonds.li3 In composite materials, for small content, 
boron tends to form a linear chain, while as its con- 
tent increases it can Jarm structures ranging from two- 
to three-dimensional £3 As far as isolated clusters con- 
taining boron are concerned, experimentally, they have 
been studied_.by Anderson and coworkersJij and by La 
Placa et alo In addition, Chopra et al£3 and Lee et 
al.Ej have experimentally synthesized the boron-nitrogen 
nanotubes and cage-like boron-nitrogen structures. Sev- 
eral authors have also i^papLciL ^carj^i^aLcalrul atopa 
on the boron ^usters.yyia'ByMBBBO'BBB'B 
Boustani et al.El have shown theoretically that, similar 
to carbon, boron has a strong potential to form stable 
nanotubular structures, Boustani and coworkers recently 
studied small cationiccHI and neutral boron clusters^ and 
obtained structures that are fundamentally different from 
crystal subunits of the well-known a— and (3— rhom- 
bohedral phases of boron, which consist mainly of B\2 
icosahedra. They classified the boron clusters into four 
topological groups: cojayex and sphericalEj , quasipla- 
narc3 and nanotubularEl The quasiplanar and convex 
structures can be considered as fragments of planar sur- 
faces, and as segments of hollow spheres, respectively. 
The main focus of their theoretical work has been to as- 
certain the structures of larger boron, c l u sters in terms 
of a small number of building blocks .EZrc3 However, re- 
cently Sabra and BoustaniE3 studied the ground-state en- 
ergetics of ladder-like quasi-one-dimensional clusters of 
boron by quantum chemical methods. They concluded 
that such structures are not the lowest in energy. How- 
ever, because of the proximity of their energy to that of 
the true ground-state-geometries, they can be regarded 
as metastable states.EH Thus, with some experimental 
manipulation, it may be possible to realize such struc- 
tures in laboratory. Keeping this possibility in mind, we 



decided to compute the static polarizabilities of ladder- 
like quasi-one-dimensional structures of boron by ab ini- 
tio many-body methods. In addition to the quantum- 
chemical calculations, we also perform the DFT-based 
calculations of static dipole polarizabilities of these clus- 
ters using the same basis set, so as to understand the 
influence of electron-correlation effects on the polariz- 
abilities of these systems. Recently, Reis and cowork- 
ers computed the static dipole polarizabilities of_Lhom- 
bic B4 using various quantum-chemical methods ,E3 and 
several other boron clusters B„(n = 3 — 8,10) within 
the framework of DFT-cmploying a variety of exchange- 
correlation functionals £9 However, unlike the quasi-orve-. 
dimensional geometries considered by us, Reis et al.E2l 
performed these calculations on the ground-state geomfe 
tries of Boron clusters optimized earlier by Boustani.EZI 
We compare our many-body static polarizabilities of-s 
ious boron clusters to those reported by Reis et alB3 
in order to understand the influence of the cluster struc- 
tures on their static polarizabilities. 

The remainder of the paper is organized as follows. In 
section [n], the applied methods and computational de- 
tails are briefly described. The results are then presented 



presented in section IV 



and discussed in section HI. Finally, our conclusions are 
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METHODS AND COMPUTATIONAL 
DETAILS 



For the closed-shell clusters, first the polarizabili- 
ties are calculated by using the restricted Hartree-Fock 
(RHF) method, and thereafter the electron-correlation 
effects are included via the M0ller-Plesset second-order 
perturbation theory (MP2), and the coupled-cluster 
(CC) techniques. For the open-shell clusters, the calcula- 
tions are initiated by the multi-reference self-consistent - 
field (MCSCF) method, while the electron-correlation 
effects are taken into account by the multi-reference 
configuration-interaction (MRCI) method. In order to 
calculate static dipole polarizability first we performed 
calculations without an external electric field, and then 
we added an external electric field of strength 0.001 a.u. 
along the x, y and z axis separately. Stability of the 
results with respect to the value of the field was care- 
fully examined by performing some calculations for var- 
ious other values of the field strength. However, when 
we perform high-level correlated calculations, the expec- 
tation value of the dipole moment is not directly avail- 
able. Therefore, to calculate the static dipole polarizabil- 
ities, we have adopted a finite-difference formula in which 
the diagonal polarizability tensor elements are obtained 
through the second derivative of the total energy with re- 
spect to the external electric field. The field-dependent 
total energy is used in the following finite-difference for- 
mula: 
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where E{ej) is total energy with respect to field £j — 

0. 001 a.u. and E(0) is total energy without field. This 
equation holds only for centrosymmetric systems. 

For the linear carbon chain, assuming the z axis is the 
chain direction, we calculated the parallel (a zz ) and per- 
pendicular (a xx ) components of the static dipolc polariz- 
ability. Our calculations arenperformed using the geome- 
try reported by Watts et alia. Since the ground state of 
even number of carbon atoms is triplet and of odd num- 
ber of carbon atoms is singlet we calculated the static 
dipole polarizabilities of even number of carbon atoms 

1. e., O4, Cq, C$ and C%o for its ground state by the MC- 
SCF and the MRCI methods, whereas for odd number of 
carbon atoms i.e., C3, C5, C7 and Cg we calculated them 
by using the RHF, MP2, CCSD, and the CCSD(T) meth- 
ods. All calculations were performed with the MOLPRO 
molecular orbital ab-initio program packageEH by employ- 
ing Sadlej basis setsE3 which was specially constructed for 
the calculation of dipole polarizabilities. 

For the chain-like boron clusters we assumed that the 
boron atoms were lying in the xy plane, with the chain 
direction along the x-axis. We first optimized the geom- 
etry of each cluster, i.e., B 4 , B e , i? 8 , Bio, B i2 and B14 
for its ground state at the B3LYP/6-31+G(d) lev&Lof ap- 
proximation by using the GAUSSIAN98 programed. The 
ground state is singlet for B4, B w , B 12 -B14 and triplet for 
£?6, B&. Then we calculated the parallel (a xx ), transverse 
(a yy ) and perpendicular (a zz ) components of the static 
dipole polarizabilities with standard polarized valence 
double-zeta (VDZ) basis sets at the Hartree-Fock and 
correlated level e.g., MRCI, MP2, CCSD and CCSD(T) 
by employing .MOLPRO molecular orbital ab initio pro- 
gram packaged. Although the basis set which we used 
in chain-like boron clusters is a rather small basis set, a 
larger set would have been computationally too expensive 
when we prolong the chain. It has been shown further 
in Ref.EII that using the larger triple-zeta basis set aug- 
cc-pVTZ does not have a large effect on the calculated 
polarizabilities for the cluster B 4 . Therefore, for these 
clusters the chosen basis set should be sufficient. 



in the same figure we. have also plotted the DFT-based re- 
sults of FucntealbaE!]. It is not possible for us to compare 
oijix re^iBw)bcreS|)^)ments because of the absence of any 
pgiacizabilityjliata on the carbon chains. From Fig. |^ it 
is obvious that parallel component of the static dipole 
polarizability per atom, a zz /N, increases roughly lin- 
early from C2 to C10, whereas perpendicular components 
ctxx/N and a yy /N are essentially constant as a function 
of N. As far as the comparison of our results with the 
DFT results of FuentealbaHl is concerned, the agreement 
is generally very good on all components of polarizabil- 
ities except for the case of N = 2. For N = 2, how- 
ever, FuentealbaE!! reports an anomalously large value 
of a xx (and hence a xx /N), making it in disagreement 
with his values of a xx /N computed for higher values of 
N . From Fig. [2] it is also clear that the directionally- 
averaged polarizability per atom a av /N, also shows an 
overall increase as a function of the chain length. Since 
polarizability is an extensive quantity, therefore, for very 
large number of atoms in the chain (N — » 00), a av /N 
should approach its bulk value. However, from our re- 
sults it is obvious that for N = 10, a av /N is still in- 
creasing as a function of N, exclusively because of the 
increase in the parallel component a xx /N. One can un- 
derstand the increase in a xx /N as a function of N on 
the intuitive grounds based upon the behavior of 7r elec- 
trons. 7r electrons (of which the carbon chain has two per 
atom), as against the a electrons, are highly delocalized 
along the chain direction. Therefore, it will take much 
larger cluster sizes before their response to an external 
field approaches that in the bulk. 

The other somewhat surprising aspect of our results 
for the carbon chains is the generally excellent agreement 
obtained between the DFT values, and the many-body 
values of the static polarizabilities. This means that for 
the static polarizabilities of carbon chains, DFT is able 
to describe the electron correlation effects quite well. It is 
also rather interesting to note that MP2 method provides 
a theoretical description of these clusters quite close to 
that obtained with the CCSDIT) method. A similar ef- 
fect was observed by Marouliso in the polarizability cal- 
culations of a system composed of two water molecules. 



B. Chain-like boron clusters 



III. RESULTS AND DISCUSSIONS 

A. Linear carbon clusters 

The calculated Cartesian components of static dipole 
polarizabilities and average polarizabilities a av = (a zz + 
2&xx)/3 for linear carbon clusters are presented in Table 
H Figure || presents our calculated polarizabilities per 
atom, plotted as a function of the number of atoms in 
the chain (N). Additionally, for the sake of comparison, 



Earlier Sabra et al.c3 had shown that strictly one- 
dimensional chains of boron are unstable. They demon- 
strated that boron prefers to form—a zig-zag ladder- 
like quasi-one-dimensional structure.Ej Therefore, in the 
present work we have concentrated on the identical struc- 
tures of boron B2„ (n = 2 — 7), which, as shown in Fig. 
|l|, can be obtained by adding boron dimers to B4, which 
has a parallelogram structure. First we optimized the 
ground-state geometry of each of these clusters by em- 
ploying B3LYP/6-31+G(d) method, and these optimized 
geometries are given in Fig. |l|. From optimization we 
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found that the system is completely centrosymmetric. 
The optimized geometries of each cluster are comparable 
to those obtained by BoustaniE3 who optimized the struc- 
tures of elemental, convex and quasiplanar boron clusters 
B n (n = 2 - 14) at the RHF level with the 3 - 21G basis 
set. 

Results of our calculation on longitudinal (a xx ), trans- 
verse (cx V y), perpendicular (a zz ), and directionally- 
averaged polarizabilities a av = {ot xx + otyy + oi zz )/2> at 
the HF and correlated level, as well as at the DFT level 
are presented in table |l[ The polarizabilities per atom 
based upon this data are plotted as a function of the num- 
ber of atoms in the cluster in Fig. [| An inspection of the 
table and the figure reveals the following trends: (i) The 
longitudinal static dipole polarizability per atom a xx /N 
increases almost linearly with N, the component a yy /N 
shows a gradual decrease, while a zz /N exhibits satura- 
tion, (ii) For the triplet ground state clusters Bg and 
Bg, the polarizabilities computed by the MCSCF and the 
MRCI methods are in almost complete agreement indi- 
cating that the MCSCF method has already captured the 
most important correlation effects, (iii) For the remain- 
ing clusters whose ground state is singlet, the inclusion 
of electron correlation effects leads to a reduction of the 
a xx component, while the other components are rather 
unaffected. For example, the CCSD(T) value of a av for 
B14 is about 6 % smaller compared to its RHF value, (iv) 
Similar to the case of carbon chains, for all boron clusters 
considered here, generally there is very good agreement 
between the polarizabilities computed by the best wave- 
function methods (MRCI and CCSD(T)), and the ones 
computed using the DFT/B3LYP approach. Thus, in 
this case also, the DFT is able to account for the elec- 
tron correlation effects quite well. 

Since there are no experimental results for the static 
dipole polarizabilities for the ladder-like structures of 
boron, we compare our respite to the theoretical results 
obtained by other authorsBj'E2l First considering the case 
of rhombic B4, for the average static dipole polarizability 
a av , we obtained 51.99,-a.u. for with the CCSD(T) ap- 
proach while Reis et al.EJ report a CCSD(T) value 60.00 
a.u. for the same quantity. The a xx and a yy values from 
our calculations cannot be directly compared to those 
reported by Reis et al.Ea because of the different orienta- 
tions of the x and the y axes in their calculations. How- 
ever ,rfor a zz , whose values can be compared directly, Reis 
et al£3 report the value 39.5 a.u., while we obtained 27.12 
a.u. for the same quantity. Although the optimized ge- 
ometries, as well as the basis set used by Reis et al.ca were 
somewhat different from ours, we still believe that those 
factors cannot explain the difference of « 12 a.u. in the 
values of a zz . However, clearly it is this disagreement — 
the reasons behind which are not clear to us — which is 
primarily responsible for the disagreement in the values 
of flu,, observed between our results and those of Reis et 
alf 

Besides B4, there are no theoretical results on the 
larger ladder-like clusters of boron. However, in another 



paper Reis et alB9 reported DFT based calculations of 
the static polarizabilities of convex and quasiplanar B„ 
(n = 3 — 10) clusters whose geometries were optimized 
earlier by BoustaniEU Therefore, in order to understand 
the effect of the geometric structure on the polarizabili- 
ties of boron clusters, we compare our results on Bg, Bg, 
and B+p, with the corresponding isomers studied by Reis 
et al.Ej For Bg Reis et al£9 considered a benzene-like 
hexagonal geometry with the Dih symmetry, and coip=. 
puted the value a av = 101.3 a.u. For B s also Reis et al.EHl 
considered a ring-like structure with the D-?h symmetry 
and reported a av = 114.6 a.u. Finally, for Bio they con- 
sidered a quasiplanar structure with the C2/1 symmetry 
and calculated a av — 143.7 a.u.E3 These can be com- 
pared with our DFT values of a av of the ladder- like Bg ; 
Bg, and Bio which were obtained to be 92.35 a.u., 130.65 
a.u., and 185.53 a.u., respectively. From the comparison 
it is clear that although the polarizability of the benzene- 
like Bg is larger than that of the ladder- like Bg, how- 
ever, for clusters containing larger number of atoms (Bg 
and Bio) quasi-one-dimensional ladder- like structures are 
more polarizable than the planar structures. Although, 
no theoretical results on a av for the planar structures of 
B12 and B14 are available, however, it is clear that even 
for those clusters ladder-like structures will be obtained 
to be more polarizable. The fact that for larger number 
of atoms, ladder-like boron chains will be more polariz- 
able than quasiplanar isomeric structures of boron, can 
be understood based upon intuitive arguments. As Sabra 
et al.cS showed by explicit calculations, with the increas- 
ing size, the 7r-electron population of these ladder-like 
chains increases. Since the 7r-electrons are quite delo- 
calized along the chain direction, their response to the 
electric fields directed along the chain direction will be 
quite large leading to large values of a xx obtained in our 
calculations. The fact that a xx increases quite rapidly 
with the increasing number of atoms also confirms this 
hypothesis. Although, the quasiplanar structure of boron 
also have ir electrons due to the sp 2 hybridization, how- 
ever, their response to the external field is distributed 
in two directions due to their two-dimensional charac- 
ter, leading to smaller polarizabilities. Therefore, it is 
the combined effect of tt electrons and the reduced di- 
mensionality which makes the ladder-like chains of boron 
more polarizable than its quasiplanar counterparts. 

IV. CONCLUSIONS AND FUTURE DIRECTIONS 

In conclusion, we have reported a systematic ah ini- 
tio study of the static dipole polarizability of the lin- 
ear carbon chains and the ladder-like boron chains em- 
ploying Hartree-Fock, many-body, and the DFT-based 
approaches. For closed-shell clusters the polarizabilities 
computed by the RHF method were generally within 
5% — 6% agreement with the ones computed by the 
CCSD(T) method. Similarly for the open-shell clusters 
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the MCSCF polarizabilities were found to be in very good 
agreement with the ones computed by the MRCI method. 
Additionally, DFT-based results on the polarizabilities 
were found to be in very good agreement with the ones 
computed by the many-body methods. This suggests the 
possibility that by employing computationally less expen- 
sive approaches such as RHF, MCSCF, DFT etc. one 
can perform similar calculations on even larger and more 
complex clusters, and obtain reasonable results on static 
polarizability. We believe that such a line of investigation 
should be pursued in future calculations. 

Our results demonstrate that in both the systems the 
component of the polarizability along the chain direction, 
as well as the average polarizability, increase with the 
chain size. This is fully consistent with presence of the 
delocalized 7r-electrons in these one-dimensional clusters. 
Thus both types of clusters should be useful in nonlin- 
ear optical applications as well. Of late, the molecular 
transport properties pfj carbon chains have been of much 
interest to physicists,E3 because of the presence of the de- 
localized 7r-electrons in them. However, our polarizabil- 
ity calculations point to the presence of the delocalized 
7r-electrons also in the ladder- like boron clusters, thus 
rendering them possibly useful in molecular-transport- 
based applications. First principles studies of the excited 
states, and the transport properties of such clusters will 
be the subject of future investigations. 
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FIG. 2. The static dipole polarizability per atom of linear 
carbon clusters. The line has been plotted to guide the eyes. 
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TABLE I. Linear carbon clusters: Static dipole polariz- 
abilities (in a.u.) calculated with the Sadlej basis set. 



Atoms 


ivietnoas 








v 12 




on nn 


9£ 99 


97 C\1 




i vi rivj i 


94 


9^ 94 


94 7Q 




T\T? r T /DOT VD a 

JJr 1 / JoO-Li r 


1 nn nn 


qq on 


7o nn 


f 1 

c 3 


I> 1_J T7 

KHr 


24.72 


51.16 


33.53 




MP9 


9/1 /I 


01. Oo 


00.00 




OUoJJ^ 1 J 


9/1 Q9 
z4.oz 


£9 Rf\ 


Q/i in 




T\T? r T /DOT VD a 


r)C on 


OU.oU 


OO. f U 


r\ 




9Q RF\ 


yi.oo 


^0 ^4 




MRCI 


29.71 


92.38 


50.60 




DFT/B3LYP a 


32.10 


93.70 


52.60 


C 5 


RHF 


35.42 


145.74 


72.19 




MP2 


35.38 


149.74 


73.50 




CCSD(T) 


35.50 


149.88 


73.63 




DFT/B3LYP" 


36.20 


142.90 


71.70 


c 6 


MCSCF 


40.81 


219.88 


100.50 




MRCI 


40.45 


220.40 


100.44 




DFT/B3LYP a 


42.80 


214.70 


100.00 


C 7 


RHF 


46.54 


303.50 


132.19 




MP2 


46.22 


314.80 


135.75 




CCSD(T) 


46.74 


312.30 


135.26 




DFT/B3LYP a 


47.00 


285.80 


127.00 


c 8 


MCSCF 


51.82 


464.78 


189.47 




MRCI 


51.34 


462.11 


188.26 




DFT/B3LYP a 


52.60 


403.00 


169.00 


C 9 


RHF 


57.58 


537.12 


217.43 




MP2 


57.20 


555.04 


223.15 




CCSD(T) 


57.64 


551.54 


222.27 



TABLE II. Chain-like boron clusters: Static dipole polar- 
izabilities (in a.u.) calculated with the VDZ basis set. 



Atoms 


Methods 


O-xx 


®~yy 


Qzz 




B 4 


RHF 


82.72 


65.54 


27.30 


58.52 




MP2 


66.34 


57.28 


27.30 


50.31 




CCSD(T) 


69.88 


58.98 


27.12 


51.99 




DFT/B3LYP 


74.60 


61.60 


29.20 


55.13 


B 6 


MCSCF 


145.45 


92.86 


40.78 


93.03 




MRCI 


144.64 


92.16 


40.28 


92.36 




DFT/B3LYP 


145.51 


91.15 


40.39 


92.35 


B 8 


MCSCF 


210.72 


101.93 


62.73 


125.12 




MRCI 


211.28 


101.17 


61.54 


124.66 




DFT/B3LYP 


229.05 


102.40 


60.51 


130.65 


Bio 


RHF 


361.86 


121.24 


97.00 


193.37 




MP2 


362.08 


120.86 


63.54 


182.16 




CCSD(T) 


345.70 


118.90 


79.82 


181.47 




DFT/B3LYP 


350.20 


118.00 


88.40 


185.53 


B12 


RHF 


527.32 


153.60 


77.90 


252.94 




MP2 


445.82 


129.52 


76.04 


217.12 




CCSD(T) 


469.70 


140.06 


76.58 


228.78 




DFT/B3LYP 


505.80 


142.20 


76.80 


241.60 


Bl4 


RHF 


698.90 


160.16 


115.36 


324.81 




MP2 


680.52 


152.68 


75.92 


303.04 




CCSD(T) 


664.48 


155.08 


97.66 


305.74 




DFT/B3LYP 


683.80 


155.20 


104.00 


314.33 



Cio MCSCF 62.85 777.51 301.07 

MRCI 62.29 776.39 300.32 

a Taken from RefM 







